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We generalize the factorized resummation ol multipolar waveforms introduced by Damour, Iyer 
and Nagar to spinning black holes. For a nonspinning test-particle spiraling a Kerr black hole in the 
equatorial plane, we find that factorized multipolar amplitudes which replace the residual relativistic 
amplitude ftm with its i-th root, pern = fel!^, agree quite well with the numerical amplitudes up 
to the Kerr-spin value q < 0.95 for orbital velocities v < 0.4. The numerical amplitudes are 
computed solving the Teukolsky equation with a spectral code. The agreement for prograde orbits 
and large spin values of the Kerr black hole can be further improved at high velocities by properly 
factoring out the lower-order post-Newtonian contributions in pi^. The resummation procedure 
results in a better and systematic agreement between numerical and analytical amplitudes (and 
energy fluxes) than standard Taylor-expanded post-Newtonian approximants. This is particularly 
true for higher-order modes, such as (2,1), (3,3), (3,2), and (4,4) for which less spin post-Newtonian 
terms are known. We also extend the factorized resummation of multipolar amplitudes to generic 
mass-ratio, non-precessing, spinning black holes. Lastly, in our study we employ new, recently 
computed, higher-order post-Newtonian terms in several subdominant modes, and compute explicit 
expressions for the half and one-and-half post-Newtonian contributions to the odd-parity (current) 
and even-parity (odd) multipoles, respectively. Those results can be used to build more accurate 
templates for ground-based and space-based gravitational-wave detectors. 

PACS numbers: 04.25.Nx, 04.30.Db 



I. INTRODUCTION 

An international network of kilometer-scale laser- 
interferometric gravitational-wave detectors, consisting 
of the Laser-Interferometer Gravitational-wave Observa- 
tory (LIGO) [l and Virgo Q are currently operating 
at the best sensitivity ever in the frequency range 10-10'^ 
Hz. We expect that in the next decade the Laser Interfer- 
ometer Space Antenna (LISA) Q will be also operating, 
but in the frequency range 10~''-10~^ Hz. 

Binary black holes are among the most promising 
sources for those detectors. During the last thirty years, 
the search for gravitational waves from coalescing black- 
hole binaries with LIGO, Virgo and LISA has prompted 
the development of highly-accurate, analytical template 
families to be employed in matched-filtering analysis. 
Those template families are based on the post-Newtonian 
(FN) approximation of the two-body dynamics and grav- 
itational radiation [1, IH. In FN theory the multipo- 
lar waveforms are derived as a Taylor-expansion in v/c 
{v being the binary characteristic velocity and c the 
speed of light). More recently, Damour, Iyer and Na- 
gar [s^l have proposed the so-called factorized resum- 
mation of the multipolar waveforms. In this resummation 
the Taylor-expanded multipolar waveforms computed in 
FN theory are re- written in a factorized, resummed form. 

In Refs. the factorized waveforms for a test parti- 

cle orbiting around a Schwarzschild black hole were com- 
puted, including also the case of comparable-mass non- 



spinning black holes. It was found that factorized wave- 
forms agree better with numerical (exact) results than 
Taylor-expanded waveforms. In particular, in the test- 
particle limit, Ref. @ compared the analytical factor- 
ized (/, m) modes and gravitational-wave energy flux to 
the numerical results obtained by Berti Q, solving the 
Teukolsky equation. The factorized waveforms have been 
also employed in the effective-one-body formalism and 
compared to waveforms computed in numerical-relativity 
simulations 0,1 [13. Also in this case, the agreement of 
the factorized waveforms to the numerical waveforms is 
better than the one of the Taylor-expanded waveforms, 
especially during the last stages of inspiral and plunge, 
and close to merger. 

In this paper we extend the factorized multipolar wave- 
forms to the case of a test-particle orbiting around a 
Kerr black hole on the equatorial plane. In the case of a 
test particle orbiting around a Schwarzschild black hole, 
the Taylor-expanded multipolar waveforms were derived 
through the FN order needed to compute the 5.5FN en- 
ergy flux [111 ] , although their explicit formulas were not 
available in the literature. In the case of a test particle 
orbiting around a Kerr black hole on the equatorial plane, 
spin terms in the Taylor-expanded multipolar waveforms 
were derived through the FN order needed to compute 
the 4FN energy flux [l2| , but their explicit formulas were 
not published. Motivated by this work, Tagoshi and 
Fujita [l^ have recently computed the spinning and 
nonspinning Taylor-expanded multipolar waveforms up 
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TABLE I: PN orders currently available in the multipolar waveforms Ctm (in the adiabatic approximation Cim = —m^ hem)- 
In the first two rows, we list the nonspin and spin PN orders beyond the leading order Newtonian term C22'''''' . In the next two 
rows, we list the nonspin and spin PN orders beyond the leading-order term for each mode Cg^''^''\ In the last two rows, we list 

the PN orders beyond the leading-order term for each mode c\^^^^ that are needed to compute the nonspin 5.5PN-energy-fiux 
and the spin 4PN-energy-flux. For each Cim, the two columns refer to the parity of the multipolar waveform ep = and = 1. 



to 4PN and 5.5PN order (see Tabic |I] for a summary), 
respectively. Also, recently, Fujita and Iyer [ij] have in- 
dependently computed the nonspinning Taylor-expanded 
multipolar waveforms up to 5.5PN order. 

Since, as said above, explicit expressions of the Taylor- 
expanded multipolar waveforms are not available in the 
literature, even at lower PN orders [lll,[l3 > we write those 
expressions explicitly in this paper (see Appendix |X| de- 
composing them in -2 spin- weighted spheroidal harmon- 
ics. Then, we apply the transformation from -2 spin- 
weighted spheroidal harmonics to -2 spin- weighted spher- 
ical harmonics, and build the factorized multipolar wave- 
forms decomposed in -2 spin-weighted spherical harmon- 
ics. The latter decomposition is the one commonly used 
in the fields of numerical relativity and gravitational- 
wave data analysis. Finally, we compare the factor- 
ized waveforms to numerical (exact) waveforms for a 
test particle orbiting around a Kerr black hole, on the 
equatorial plane, solving the Teukolsky equation [l5l - [T7| . 
Finally, we derive the factorized multipolar waveforms 
for spinning, nonprecessing black holes of comparable 
masses. Those factorized waveforms were recently used 
in the spinning efTective-one-body model of Ref. [l^l and 
compared to numerical-relativity simulations of spinning, 
nonprecessing equal-mass black holes from the Caltech- 
Cornell-CITA collaboration. 

This paper is organized as followed. In Sec. |lT] we 
work out the factorized waveforms decomposed in -2 spin- 
weighted spherical harmonics for a test-particle orbiting 
a Kerr black hole, on the equatorial plane. In Sec. IIII[ 
we compare the gravitational-wave energy flux and the 
(Z, to) modes of analytical factorized waveforms to numer- 
ical waveforms. The numerical results arc obtained solv- 
ing the Teukolsky equation [l5l - [l7| . In Scc.|lV]we derive 
the factorized waveforms for generic mass-ratio spinning, 
non-precessing black holes. Section |V] summarizes our 
main conclusions. In Appendix |X] we write the Taylor- 
expanded multipolar waveforms in the test-particle limit 
through the PN order currently known. In Appendices iBl 
[Cl IdI and |E] we give the complete expressions of the 
/fm's, Cf„i's, pfm's and Stm's, for 4 < / < 8. Finally, in 



Appendix |F] we compute the I and to dependence of the 
spin terms in the mass and current multipole moments 
at 0.5PN order and 1.5PN order, respectively. 



II. FACTORIZED MULTIPOLAR WAVEFORMS 
FOR A TEST PARTICLE ORBITING AROUND A 
KERR BLACK HOLE 

We consider a nonspinning test-particle orbiting 
around a Kerr black hole, and extend the factorized wave- 
forms of Ref. to the case where the motion is quasi- 
circular and confined to the equatorial plane, that is the 
spinning, nonprecessingcase. The factorized multipolar 
waveforms are built as [g 

htm = '"^ him = '"^ 'S'cff ^ Tfm e"^''" fgm ■ (1) 



where denotes the parity of the multipolar waveform. 
In the quasi-circular case, Cp is the parity of £ -I- to: = 

7r(£ + to). The leading term /i^^''"^ in Eq. ^ is the 
Newtonian order waveform 

'^^"'"^-^"te^..e.(.).(^-^)r---™ , 

(2) 

where v is the orbital velocity, y^™(6',0) are the scalar 
spherical harmonics, nf^ are 
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and ce+ep{v) are functions of the symmetric mass ratio 
v = mi TO2/A/^, with M = mi + TO2: 
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Although in this section wc consider the test particle limit 
TOi = M 3> m2 = fJ., that is — ^ 0, the above relations 
will be used for generic v in Sec. IIVI and Appendix iFl 

We shall define the source factor S'^^'''' and the tail fac- 
tors Tgm in Sec. Ill Bl In Sees. Ill Bl III CI we compute 
the imaginary and real PN spin effects in the e"'*™ 's and 
/fm's respectively. We shall obtain those quantities by re- 
quiring that when we Taylor expand the factorized wave- 
forms dH) the results coincide through 4PN order, for the 
spin terms, and 5.5PN order, for the nonspinning terms, 
with the Taylor-expanded waveforms given in Sec. Ill Al 
and Appendix 1X1 



where a is the spin of the Kerr black hole, having the 
dimension of length (while we also define q = a/M), and 
Wo = mil is a multiple of the orbital frequency f2. The 
-2 spin- weighted spherical and spheroidal harmonic bases 
are related by 



(6) 



The coefficients and are given in Ref. [i2| as 



A. Taylor-expanded multipolar waveforms 

The Newman- Penrose scalar ^'4 = — — ih^) can 
be decomposed in either -2 spin-weighted spherical har- 
monics -2Yem{0, 4>) =-2 Ptm{0) e™"^, or -2 spin- weighted 
spheroidal harmonics -2SZ°{0) e™"^ as 



2 / (£+3)(£-l)(f+m+l)(£-m+l) 



ie+2)(e~2)(e+m.){e- 
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£' = e-i 

otherwise 
(7) 
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<em\ 



e rn=-e 
= C -^Pernio) 



jLJo{r'-t) 



2 J ImyJ) iujoir' ~t)+im<j> 



/27r 



and if £' = f we have 



= J2 J2 ^fa,x.o e»c.o(r--t)-H»m0^(5) while if £V ^ we have 



i'^Im ) + i^im ) 



(8) 



-4c^+S/ + 1, 1, 0|^', m) {£ + 1, 2, 1, 0\£\ 2) - 4c 



e-i /2£-l 



- \j ^ir^i^^ ™' 2, 0|£', m) (£, 2, 2, 0|£', 2) 



l,TO,l,0|f,m)(£- 1,2, 1,0|€',2) 

(9) 



where (ji, mi, ja, ^2!^, M) is a Clebsch-Gordan coeffi- 2 (£ 4)(£ - 3)(^^ £ - 3to^) 

cient and 3 1)(2£ -f 3)(2^ - 1) ' ^ ' 

Xo{£) = + 2), (10) 

Ai(^,m) = -2m(^2+^ + 4)/(^2_^£), (11) . . 

2 2 2 nonspmmng case, _20^„ (p) reduces to 

A2(^, m) = -2(£ + 1) (c^,t/) + 2£ (c^~^) + - -2Pem{0) which has the closed expression 
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The -2 spin-weighted spherical harmonic basis is or- 
thonormal in the sense that 




-2Pe,niO)-2 Pi'm' (9) sin dO = See, Smm' ■ (14) 



To exphcitly write the modes Cem and ZimujQ expanded 
in w, we find it convenient to introduce the following no- 
tation 

ZlmuJo = ^^'to^q' ^fmwo > (16) 

where C^^''^"^ and zj,^'^^^ represent the Newtonian con- 
tributions and, as said above, ep denotes the parity of 
the multipolar waveform. In the adiabatic limit, C^m = 
—m^n^him. Therefore, whereas the Newtonian contri- 
bution to Cfm and hem differ by a factor of —m^fl^, 
the PN corrections are the same, i.e., Cein = hem- The 
Newtonian contributions in the Cem 's or Zfmuia are [see 
Eq. @], 



xr^-<^'"-"(7r/2,0), (17) 

where we define i; = {MVLfl'^. In Refs. [H d, the 
Taylor-expanded multipolar waveforms were calculated 
at the PN order needed to compute the nonspinning 
5.5PN-energy-flux and spin 4PN-energy-flux, respec- 
tively. For the purpose of the present paper, Tagoshi 
and Fujita [ll] extended the computation of the mul- 
tipolar waveforms at higher PN order. Although those 
new PN corrections are not sufficient for computing the 
energy flux at the next order (6PN and 4.5PN order in 
the nonspinning and spinning cases, respectively), they 
do improve our knowledge of the multipolar waveforms, 
as we shall discuss below. In Table U we list the PN 
orders available to us in each multipolar waveform C^m, 
while the explicit Taylor-expanded waveforms, Zemu)a^i 
arc given in Appendix [XI 

We compute the Cem^s from the Zemujo^^ by applying 
Eq. (jni) and the orthogonality condition of the -2 spin- 
weighted spherical harmonics 
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+ aUJO ^ 4'rn -2-P^"m + (a Wo)^ ^ 4'm 



ZemujO + a^^a C^e'm Zl'muO + (a^^o)^ '^('m Ze'mujQ + C( 



(18) 



We notice that the mixing of spheroidal waveforms hap- 
pens among modes with the same m and different £. 

The Cem raodes are computed in perturbation the- 
ory [nl, [H, [1^, [iOl using a coordinate swtem different 
from the one used in PN calculations [21I, When 
expressing both modes in terms of the orbital frequency 
they should coincide. However, the presence of tail terms 
in both calculations demands a careful treatment. In PN 
calculations, the tail terms contain a freely specifiable 
constant rp that corresponds to the difference in the ori- 
gins of the retarded time in radiative coordinates and in 
harmonic coordinates in which the equations of motion 
are given (see e.g., Eq. (3.16) in Ref. [22|). This constant 



can be absorbed into the phase of the PN modes (see 
e.g., Eq. (8.8) in Ref. [22|) once it is traded with xq (or 
vo) mi as 

11 2 4 2 

log EE 2 log Wo = — - - jE - - log 2 - - log(ro/M) , 

(19) 

where 7^; = 0.577 215... is the Euler's constant, and 
throughout the paper, we use "log" to denote the nat- 
ural logarithm. In perturbation theory calculations, 
Schwarzschild or Boycr-Lindquist coordinates are used. 
The waveforms at infinity are naturally expressed with 
the tortoise coordinate, and the relation between the 
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Schwarzshild or Boyer-Lindquist coordinate and the tor- 
toise coordinate has an arbitrary constant which in 
Refs. [lll,[il,[ll[23 is fixed to -2M log(2M). [H 

We find that to recover the PN results we need to 
express some of the 7£;'s in the perturbation-theory 
modes [31 in terms of xq and ro using Eq. ([T9| and set 
To = 2M/y/e. More specifically, we replace some of the 
7£;'s using the following equation [2l| 

log«o = ]^~]-iE-\\og2-\ log(2e-i/2) . (20) 
36 3 3 3 

We notice that the constant 7'o will appear later in our 
definition of the tail term T^m of the factorized resummed 



waveforms. In fact, since the Tirn term resums all tail in- 
tegrals that contain rg at known orders, it is the only 
term in the resummed waveforms that depends on tq. 
Finally, to ease the notation, we follow Ref. and in- 
troduce eulerlog„(w^) = 7e -I- log2 -I- logm -I- l/21og?;^ 
into our Cgrn expressions. 



Below we list the Cgrn^ through I = 4, and give the 
expressions for 4 < £ < 8 in Appendix |B] The differences 
between the Cim's and Z^rnuio's concern only spin terms. 
We obtain 
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where the Z^j^s can be found in Appendix 1X1 We no- (relative to Z^^s), the Cim's do not, except for C2i- 
tice that whereas the Zim^s contain 0.5PN spin terms The spin terms in the multipolar waveforms (j21aP " (|21ip 
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agree with the currently known spin terms computed in 
Ref. [11], and with the 0.5PN and 1.5PN spin terms in 
the odd and even-parity modes computed in AppendixlF] 



B. Source and tail terms 

In the Umit of a nonspinning test particle of mass /i 
orbiting around a Kerr black hole of mass M in a quasi- 
circular equatorial orbit, the energy and orbital angular 
momentum, in Boycr-Lindquist coordinates, read [2J] 



Eir) _ l-2Af/r + aMi/Vr3/2 

^ v/1 - Ml/r + 2a Af 1/2/^3/^ ' 
L{r) _ pr l-2aA/VVr3/2 + a2/r^ 

Tm ~ V M yr^7\//r + 2a M^^T^ 



(22) 
(23) 



where r = (1 — av"^)^^^ /v"^ . The source term in the fac- 
torized waveform ([T]) is 



S. 



off 
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M/v) 



(24) 



where ^M/v is the Newtonian angular momentum. We 
use the resummed tail factor Tgm given in Eq. (19) of 
Ref. 
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where k = mCl, k = H^caik, and the real Hamiltonian 
in the test-particle limit reduces to iJrcai = ^I- Once 
again, wc emphasize that the constant vq must take a 
fixed numerical value, 2Mj ^ [s^, to reproduce the cor- 
rect test-particle limit waveforms. We notice that there 
is no spin contribution to T^m since the latter resums 
the corrections to the waveform when traveling through 
a long-range Coulomb-type potential generated by the 
mass M [25l, [2^ . Spin effects generate a short-range po- 
tential, thus they do not contribute to Tfm- 

Wc compute the phase correction factors e'**™ in 
Eq. ([1]) by Taylor expanding the factorized waveforms 
him given in Eq. ([1} , comparing the result with the Cira 
waveforms derived in Sec. Ill Al fin the circular-orbit, adia- 
batic approximation Ctm ~ —(m£t)^ him)-, and collecting 
all imaginary terms into Sim- We obtain 
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Notice that the nonspinning terms in the Sim already 
appeared in Ref. @, except for the terms at 3PN order 
(fS) [Ref. [il did compute S22 at 3PN order]. Wc find that 
those 3PN-order terms in the Sim are necessary to obtain 
full agreement between the factorized waveforms and the 
nonspinning Cim waveforms through 3PN order. Wc note 
that the nonspinning terms at 3PN order in the Sim's are 
the same as the 3PN phase terms in Zimujo [3- 
This happens because in the test-particle limit the PN 
expansion of Tim does not contain imaginary terms at 
3PN. Thus, for q = 0, the phase corrections Sim at 3PN 
order do not contain any additional terms other than the 
3PN phase terms in Zimua- We further note that some 
of the above S's can be obtained directly in the standard 
PN and test-particle limit calculations. For example, the 
terms proportional to ttv^ and tt^w^ (for q = 0) are the 
same as the phase factors in the asymptotic amplitude in 
the test-particle limit calculations (e.g., Eqs. (30)-(32) in 
Ref. [H and Eq. (4.17) in Ref. Hj). 



C. Taylor-expanded residual terms and their 
resummation 

In the circular-orbit, adiabatic approximation Cim = 
— {mfl)^ him- By Taylor expanding the factorized wave- 
forms him given in Eq. ([1]), comparing the result with the 
Cim waveforms derived in Sec. Ill Al and factoring out the 
imaginary terms in the Sim of Eqs. (|26ap - (|26ip . we derive 
the fim's in Eq. ([T]). We notice that in the case of even- 
parity modes, the determination of the fim is unique. In 
the case of odd-parity modes, it depends on the choice of 
the source which, as explained above, can be cither the 
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energy or the angular momentum. We denote with 
and the odd-parity modes computed with the energy 
and angular-momentum sources, respectively. [Since in 



both cases the source is a real quantity, the phases Sim's 
remain the same.] We obtain through £ = 4 (see the 
Appendix [Cl for modes with 4 < ^ < 8) 
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where, as introduced above, we have defined 
eulerlog„('(;2) = log 2 -I- logm -I- 1/2 log i;^ with 

m = 1,2,3,.... Note that all the nonspinning terms in 



Eqs. p7ap ~ p7i|) appear at even powers of v, and the 
spin terms at odd powers of v. Moreover, except for the 
(2,1) odd-parity mode, all the other odd-parity modes 
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do not have a spin contribution at 0.5PN order. This is 
consistent with the results of Appendix [Fl 

As emphasized in Ref. @, the decomposition of the 
Taylor-expanded multipolar waveform into several fac- 
tors [see Eq. ([1])] is in itself a resummation procedure. 
In fact, the factorization of T^m has absorbed powers 
of mTT, which introduce large coefficients in the Taylor- 
expanded waveform. Moreover, in the quasi-circular 
case assumed here, the factorization of the energy or 
angular-momentum sources, has extracted the pole lo- 
cated at the light-ring position v = ^fMjri^ with ~ 
2M [1 + cos[2/3 arccos(=Fa/-M")]] (where =F refers to pro- 
grade and retrograde orbits, respectively) which causes 
the coefficient of u^" in any PN-expanded quantity to 
grow as rj" as n — >■ oo. As we shall see in Sec. IIIH de- 
spite those improvements, the /fm's above are not close 
enough to the exact results for large velocities. 



As wc shall discuss in detail in Sec. IIIIl the /f„i's can 
be further improved by applying the p-resummation and 
Fade-summation @. The former consists in finding the 
new functions pim such that 

fim = (pimY ■ (28) 

The motivation for the p-resummation is to reduce the 
magnitude of the IPN-order nonspinning coefficient in 
fim, which grows linearly with £ (see Sec. IID of Ref. [6|]). 
We find that such dependence on £ does also affects the 
1.5PN spin terms in the even-parity modes computed as 
function of £ and m in Appendix iFl In fact, we find that 
/i^^n ^ -2£qv^/3, and so /|™" = -2£qv^/3. Thus, we 
apply the p-resummation also to the spin terms, and find 
(see Appendix [D] for modes with 4 < ^ < 8) 
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FIG. 2: We plot the pimS extracted from the numerical data as function x = v'^ . The upper panels (blue colors) refer to 
q = 0.95, the lower panels (red colors) to q = —0.95. The variable x ranges between < a; < a;LSo(cs). 
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Lastly, we may use ^(r) instead of |i| as the source term 
in Eq. ([1]) for the odd-parity modes. The corresponding 
firn ^^'^ pfm Gxpressious are given in Appendices [C] and 
IdI respectively. 

In the next section we shaU investigate the numeri- 
cal (exact) pim^s, and compare them with the analytical 
ones. We shall find that the agreement between the nu- 
merical and analytical pim is quite good, except for some 
modes. Guided by the comparison with numerical re- 
sults, we shall apply the Padc summation to the pim^s, 
and also work out an improved resummation which con- 
sists in factoring out the lower-order PN terms in the 
Taylor-expanded pirnS. We find that this factorization 
brings the zeros of the analytical pim closer to the nu- 
merical (exact) ones. 



III. COMPARISON BETWEEN ANALYTICAL 
AND NUMERICAL RESULTS FOR THE 
TEST-PARTICLE LIMIT CASE 

We have two goals to achieve in this section. The first 
is to accurately model the amplitude of the (Z, m) modes 
for several values of the spin parameter q and velocity 
V. The second is to obtain the best agreement between 
the numerical (exact) and analytical energy fluxes with- 
out introducing adjustable parameters in the analytical 
model. 

The numerical values of the energy flux used in this 
paper are obtained with a high precision numerical code 
which solves the Tcukolsky equation [l5l - [T7| . The ho- 
mogeneous solution of the radial Tcukolsky equation is 
obtained numerically by usirig a formalism developed by 
Mano, Suzuki and Takasugi [2^. In this method, the ho- 
mogeneous solutions are expressed in terms of series of 
two kinds of special functions, hypergeometric functions 
and confluent hypergeometric functions. In Refs. [islflGj . 
it was shown that the series converges very fast and one 
can compute numerically the homogeneous solutions very 
accurately. The homogeneous solution obtained with this 
method was applied to the numerical calculation of grav- 
itational waves emitted by a particle in a quasi-circular 
and equatorial orbit around a Kerr black hole [l^, [l^ . 
In this paper, for the comparison with analytical formu- 
las, we compute the Zimujo (and thus the Cim) as well as 
(dE/dt)ern for various q and O. The computation is done 
with the double precision accuracy, and the estimated ac- 
curacy of Zimuja (and thus the Ctm) as well as {dE/dt)im 



is about 14 significant figures. As in Ref. [l^, the accu- 
racy is estimated by comparing the energy fiux with that 
of Ref. [2^ in which the accuracy was estimated as about 
20 significant figures. 

A. Hierarchy between the (Z,m)'s modes 

In Fig. [T] we study the hierarchy among the 
numerically-computed modes and plot |/i^m|/|^22| ver- 
sus V for the representative spin cases: q = 0.95, 0, and 
—0.95. The parameter v varies between 0.1 and ulso('?)- 
where we denote with LSO the last stable orbit for a 
test-particle in the Kerr geometry. 

The strain waveforms h^mS are computed from the 
Cfm's under the quasi-circular adiabatic assumption, i.e., 
him = —Cgra/i'm^y'- As wc shall discuss in Sec. IIII Cl 
the energy flux for quasi-circular adiabatic orbits can be 
computed through the well-known relation 

^(^) ^ i E E (™")' ■ (30) 

£ m=-(. 

Thus, when analyzing the contribution of the hi,nS to 
the energy flux, we need to remember that = 
iCirn/imD,). Thus, the dependence of /i^m's on m is 
different than the one of /i^m's, and, as a consequence, 
the hierarchy of the modes in the energy flux is different. 

Denoting by |/ifm 1/1^22 1 the relative strain amplitude 
and by |ft.£mP/|/j22p the relative radiation power, we flnd 
the following trends. In the anti-aligned case q = —0.95 
and the nonspinning case, the (3,3), (2,1) and (4,4) 
modes are the largest subdominant modes in terms of 
strain amplitude. In terms of radiation power, they 
are also among the largest subdominant modes although 
their hierarchy changes. The (4, 4) mode contributes 
more power than the (2, 1) mode because of its larger 
TO. For the same reason, in the nonspinning case, the 
(5, 5) mode contributes more power than the (2, 1) mode 
and becomes the third strongest subdominant mode. In 
fact, in the anti-aligned and nonspinning cases, relative 
to the (2, 2) mode, the (3, 3) mode contributes > 10% 
of radiation power at the LSO, only the (3,3) and (4,4) 
modes contribute > 1%, and the (5,5) mode contributes 
1% only in the nonspinning case. In the aligned case 
q ~ 0.95, we plot in Fig. [1] the relative strain ampli- 
tudes of 8 modes that are larger than 5% at the LSO. 
In terms of the relative radiation power, the (3, 3), (4, 4), 
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FIG. 3: We plot the pi^s extracted from the numerical data as function of a; = v'^ . The upper panels (blue colors) refer to 
q = 0.95, the lower panels (red colors) to g = —0.95. The variable x ranges between < a; < Xhso{a)- 



(5,5), (6,6), (7,7) and (8,8) modes are the largest sub- 
dommant modes. The (3, 3), (4, 4) and (5, 5) modes each 
contributes > 10% relative to the (2, 2) mode at the LSO. 
In particular, the (3, 3) mode contributes > 30% relative 
to the (2, 2) mode to both the strain amplitude and the 
radiation power. Accurate modeling of its amplitude is 
therefore crucial in modeling the full gravitational-wave 
waveform and the energy flux. 



B. Comparison between the analytical and 
numerical modes 



We now examine the amplitude agreement of the nu- 
merical and analytical waveforms, focusing mainly on the 
dominant modes: (2,2), (2,1), (3,3), (3,2), (4,4) and 
(5,5). 

In Figs. [5] and [3] we show several numerical pim^s ver- 
sus X = v'^ for three representative spin cases: q = 
—0.95,0,0.95. Since the latter are real, the numerical 
pern's are obtained using Eq. ([1]) with fgm = pim^ that is 
dividing the numerical by i\Tem\ S^s)^^^ ■ The 

numerical him sltc computed from the numerical Cem 
through the relation hem = — C'fm/(wf2)^ 

Using the 0.5PN (1.5PN) order spin terms in the odd 
(even)-parity modes computed in Appendix |F] for generic 
i and m, and the non-spinning IPN terms derived in 
Refs. [i,[2l|, we have 



(31) 



rcven 



{£ + 9) 



2£ {£+l){2i + 3) 



and 



1 - - (7 2-^/2(5^2 <5,nl 



1 2 

1 + 



m2 {I + 4) 



+0(x3/2) . (32) 

Note that the 1.5PN spin terms in the odd-parity modes 
are not known for generic i and m, but they are available 
through £ — 6 in this paper. 

Reference @ pointed out that because the IPN order 
term in the /|™" and scale as £ and is negative, for 
large £ it can cause the fgm to go to zero even before 
reaching the LSO. For example if we consider the LSO 
in Schwarzschild, a;Lso(0) = 1/6 (ulso(O) = — 
0.4082), fee at IPN order has a zero at t; = 0.3634 Q. In 
the even-parity case, the inclusion of the 1.5PN spin term 
with q > can cause the zero to occur even at smaller 
values of v. In particular, for q ~ 0.95, /ge has a zero at 

V = 0.3362 (t;LSo(0.95) = 0.6497). By contrast, the cases 
with g < can push the zero to negative or imaginary 
values, or to values of v above the LSO, thus making it 
harmless. For example, for q ~ —0.95, fee has a zero 
at w = 0.4075 (wLSo(-0.95) = 0.3373). Similarly, when 
considering the odd-parity modes for large £, e.g., the 
/(^ mode, we find that in the non-spinning case the IPN 
term causes f^^ to have a zero a.t v = 0.3602, and the 
inclusion of 1.5PN spin term causes the zero to move to 

V = 0.3502 for q = 0.95, and to v = 0.3717 for q = -0.95. 



12 



In the spin case, the above problem can be even worst 
than in the non-spinning case for lower values of I. For 
example, the IPN term causes a zero in the /as at w = 
0.5345 which is above wlso(O), but the inclusion of the 
1.5PN spin term moves the zero to w = 0.4764 for q = 
0.95 which is quite below wlso(0.95). 

Motivated by the above discussion and the result in 
Appendix [F] that shows that the even-parity 1.5PN spin 
terms scale as £ ifim'^ ~ —2lqv^/3), we adopt the p- 
resummation also for the spin terms. The pirnS through 
1.5PN order read: 

- i ^ + 2£(£+l)(2£ + 3)y' 

-\qx^''- + 0{x^) , (33) 



and 



2 I 



1 2 m^il + A) 

1 + T - ^ + ^TTT-. 1 I X 



£ £2 2^(^ + 2) (2^-1-3)^ 
. (34) 

We notice that the IPN and 1.5PN terms in pge cause a 
zero at w = 0.8902 for q = 0, and at w = 0.7577 for q ~ 
0.95. The zero in P33 occurs at v = 0.9258 for q = 0, and 
at w = 0.7765 for q = 0.95. All these numbers are larger 
than vi^soil)- Note however that the p-rcsummation may 
be less effective for q > 0.95, since at g = 1, the zero in 
P66 occurs at v = 0.7530 and the zero in P33 occurs at 
V = 0.7713, both smaller than z;lso(1) = 0.7937. Of 
course all this discussion does not take into account the 
higher-order PN terms which can also move the zero to 
lower or higher values. However, as we shall see below, 
the behavior of the numerical pim is captured by the 
0.5PN, IPN and 1.5PN terms. 

In Figs. Hand El we plot the £ = 2,3,4,5,6 (m = 
£,£ — 1, . . . , 1) numerical modes versus x. First, as ob- 



served in Ref. [6| for the nonspinning case, also for the 
spin case, the behavior of the pem is reasonably sim- 
ple. In particular, except for the (2, 1) case which shows 
a special shape due to the presence of the 0.5PN term 
i^/x), all the other modes are well represented by (bro- 
ken) straight lines with one or two changes in the slope at 
high frequency. As in the nonspinning case, but less pro- 
nounced here, for each value of £, the (negative) slopes 
of the dominant m — £ (even-parity), and subdominant 
m = £—1 (odd-parity) modes are close to each other, and 
these slopes become somewhat closer as £ increases. This 
property is reproduced by the analytical pim^s truncated 
at 1.5PN order through £ = 6 modes, whose 1.5PN terms 
are known. 

In Figs. |4] and [5] we compare the numerical and analyt- 
ical P22 and P33, respectively. We use the following nota- 
tion for the analytical models. We indicate with TN[/C^m] 
the Pim expanded in Taylor scries of v through . We 
indicate with PN[P£m] the Padc-summcd pfm with M 
and N denoting the order of the polynomial in v in the 
numerator and denominator, respectively. When apply- 
ing the Pade summation in presence of logarithms (i.e., 
log(i;)) we treat the latter as constants. We indicate with 
an improved resummation of the Taylor-expanded 
Pern's which consists in factoring out their 0.5PN, IPN 
and 1.5PN order terms, that is we write 



J 



where the coefhcients c 



Jrn 



1/2' 



r,im „ J Jm 
-1 '^'^^ ^3/2 



4"'v' + ■■■), 

(35) 

are the 0.5PN, 
IPN and 1.5PN order terms in the p^m, and the coeffi- 
cients df™ with i > 2 in Eq. ((35l) are obtained by im- 
posing that the Taylor-expanded coincides with pim- 
We shall motivate the introduction of the p£„'s in the dis- 
cussion below, but basically we find that the first factor 
on the right-hand side of Eq. (|35p can capture reasonably 
well the zeros of the numerical (exact) pfm's. 

For the modes ^ < 4, we find the following p^^'s: 
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FIG. 5: Numerical and analytical psa's as functions of x = . The three panels are for spin values q = 0.95,0 and —0.95. 
The notation of the analytical pss models follows the definition in Sec. IIIIBI In the antialigned q — —0.95 case, the numerical, 
TsIpss] and pgg lines overlap, while in the nonspinning q = case, they also overlap with P* [pss]. 
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Wc notice that for a few modes, it is convenient to factor 
out even the 2PN order term. The procedure of factoring 
out zeros of pim can be improved in the future by intro- 
ducing appropriate adjustable parameters and cahbrate 
them to the numerical result. 

In Figs, m and [5] we also show results when adopting the 
Fade summation. We find that the diagonal and nearest- 
diagonal Fade-summation improve the agreement with 
the numerical results not only for the (3, 3) mode, but 
also for the (3, 1) and (4, 2) modes. An even better 
agreement for several modes is obtained when using the 
farthest-diagonal Fade-summation. However, this quite 
interesting result suffers by the presence of spurious poles 
appearing for several g values, and for this reason we will 
no longer discuss the Fade-summation in this paper. 



Finally, we observe that close to the LSO the even- 
parity modes agree slightly better to the numerical 
results than p^'s. Thus, we adopt in this paper the 
multipolar waveforms built with the p^. In Figs. IHl [7] 
and [5] we compare the Taylor-expanded, p-'^-resummed 
and numerical Newtonian-normalized multipolar ampli- 
tudes for the dominant modes. In general, the and 
p-resummed amplitudes agree better with the numerical 
amplitudes than Taylor-expanded amplitudes do, espe- 
cially for higher-order modes. More specifically, wc find 
that p-resummed amplitudes (not shown in Figs. IHl Eland 
E]) differ from the numerical ones by <0.6% up to w < 0.4 
for the (2,2), (2,1) and (3,2) modes and by <1.8% for the 
(3,3) and (4,4) modes. Their fractional difference grows 
up to ~ 1-10 at the LSO when g = 0.95. 
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FIG. 6: Upper panel: Comparison between the numerical and analytical Newtonian normalized |/i22| and |/i2i| modes for a 
test-particle orbiting around a Kerr black hole in the equatorial plane. For the numerical data and analytical models (Taylor- 
expanded and p^-resummed), we have nine curves corresponding to different spin values of the Kerr black hole. From top 
to bottom, the spins are q = —0.95, —0.75, —0.5, —0.25, 0, 0.25, 0.5, 0.75 and 0.95. Lower panel: relative fractional difference 
between analytical and numerical for the representative spin values q = —0.95, 0, 0.95. 




FIG. 7: Upper panel: Comparison between the numerical and analytical Newtonian normalized j/iaal and \hz2\ modes for a 
test-particle orbiting around a Kerr black hole in the equatorial plane. For the numerical data and analytical models (Taylor- 
expanded and p^-resummed), we have nine curves corresponding to different spin values of the Kerr black hole. From top 
to bottom, the spins are q = -0.95,-0.75,-0.5,-0.25,0,0.25,0.5,0.75 and 0.95. Lower panel: relative fractional difference 
between analytical and numerical \hirn\ for the representative spin values q = —0.95, 0,0.95. 
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FIG. 8: Upper panel: Comparison between the numerical and 
analytical Newtonian normalized \h44 \ mode for a test-particle 
orbiting around a Kerr black hole in the equatorial plane. For 
the numerical data and analytical models (Taylor-expanded 
and p-^-resummed), we have four curves corresponding to dif- 
ferent spin values of the Kerr black hole. From top to bottom, 
the spins are q = —0.95, —0.5, 0, 0.5 and 0.95. Lower panel: 
relative fractional difference between analytical and numeri- 
cal \hern\ for the representative spin values q = —0.95, and 
0.95. 



When applying the p-* -resummation, we find that 
the fractional amplitude difference between the nu- 
merical and analytical (2,2) amplitude at the LSO is 
16% (33%), 0.18% (0.32%) and 0.20% (0.85%) for q = 
0.95,0,-0.95, respectively. We indicated in paren- 
thesis the numbers when Taylor-expanded amplitudes 
are employed. For the (2,1), (3,3) and (4,4) modes, 
for which fewer spin PN terms are known (see Ta- 
ble m, the improvement due to the p-^-resummation is 
more striking. In fact, for the (2,1), (3,3) and (4,4) 
modes we obtain a fractional amplitude difference of 
2.4% (4.2), 0.2% (0.58%) and 0.0036% (0.15%), 7.5% (2), 
0.027% (0.55%) and 0.13% (0.2%), 16% (7.5), 1.7% (28%) 
and 0.6% (5.8%), for q = 0.95, 0, -0.95, respectively. 

We summarize the results of Figs. [6l [7] and |8] as fol- 
lows. First, we remark that the Taylor-expanded am- 
plitudes agree with the numerical ones quite well for the 
(2,2) mode where the PN expansion is known through the 
highest order today (5.5 PN for nonspinning terms and 
4PN for spin terms). Thus, for the (2,2) mode the im- 
provement due to the resummation technique is marginal. 
We expect that a similar result holds for higher modes 
when sufficient PN terms are known. Second, the factor- 
ized resummed waveforms consistently improve the am- 
plitude agreement with numerical waveforms for several 
values of q and large spanning of v. In the lower panels 



of Figs.|6l[7]and[8l we observe that the fractional ampli- 
tude difference between the numerical and p-'^-resummed 
waveforms is always smaller than the difference between 
the numerical and Taylor-expanded waveforms, except 
around the v values where the numerical and Taylor- 
expanded amplitudes coincide. For all modes [except the 
(2,2) mode] and all spin values shown in the figures, we 
find that p-'^-rcsummcd amplitudes are typically closer to 
the numerical amplitudes than Taylor-expanded are by 
an order of magnitude or more. 

Finally, for £ > 5 modes, the p-'' -resummation is not 
very successful in modeling the numerical amplitudes, 
but it is better than Taylor-expanded amplitudes. We 
know nonspinning and spin corrections only through 
2.5PN order in the (5,5) mode (see Table H, thus it is 
not surprising that we cannot model those modes very 
well. Since the contribution of the £ > 5 modes to the 
radiation power and strain amplitude is not negligible, it 
would be very useful to calculate higher order corrections 
in those modes in the future. 



C. Comparison between analytical and numerical 
energy fluxes 

Here we compare numerical and analytical Newtonian- 
normalized energy fiuxes for a test-particle orbiting a 
Kerr black hole in the equatorial plane. The fiuxes 
are computed by summing the power radiated using 
Eq. (PH)) and setting £ = 8. For a test- particle moving 
along a quasi-circular equatorial orbit, the Newtonian- 
normalized flux is F(v)/ Fj^cwtiv) where FNGwt(i') = 
32z/2x;io/5. 

We note that the dominant error source of the numeri- 
cal calculation of the total flux is the truncation at £ = 8 
of the mode summation. Let Fi=s{v) be the contribution 
from £ ^ 8 mode for F{v). The fraction, Fi=s{v) / F (v) , 
is about 10"^° around v = 0.1 and 10"^ to 10"^ around 
the LSO. 

In Fig. [9l we compare numerical and analytical 
Newtonian-normalized energy fluxes for different spin 
values of the Kerr black hole. In the left panel of 
Fig.|9]we consider two Taylor-expanded fiuxes computed 
from the Taylor-expanded /i^m's: one that truncates all 
terms beyond 5.5PN order and spin terms beyond 4PN 
order (Taylor-expanded-truncated) and one that keeps 
all higher order terms (Taylor-expanded-nontruncated). 
[The former is the Taylor-expanded flux that consistently 
includes nonspinning effects through 5.5PN order and 
spin effects through 4PN order [ll|, the latter in- 
cludes new hig her-order PN terms computed by Tagoshi 
and Fujita |13|.] 

In the left panel of Fig. HI we do not show the Taylor- 
expanded flux truncated at 4PN order [l2| since its agree- 
ment with the numerical flux is rather poor. Figures 2 
and 3 of Ref. [lOj show that in this case the Taylor- 
expanded flux starts to differ from the numerical one at 
a relatively low velocity of u = 0.2 for all spin values. By 
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FIG. 9: Upper panels: Comparison between the numerical and analytical Newtonian normalized energy flux for a test- 
particle orbiting a Kerr black hole in the equatorial plane. There are nine curves for each waveform's model and nu- 
merical data. They correspond to different spins of the Kerr black hole. From top to bottom, the spins are q — 
—0.95, —0.75, —0.5, —0.25, 0, 0.25, 0.5, 0.75, and 0.95. The Taylor-expanded-truncated flux includes test-particle spin terms 
through 4PN order [l^], and nonspinning terms through 5.5PN order [ll|]. The Taylor-expanded-nontruncated flux includes 
higher-order PN terms originated by the new PN terms in the hem's computed by Tagoshi and Fujita [T^] (see Table |l]). The 
p-resummed flux is plotted for one spin value q — 0.95 in the upper right panel, to show the large improvement when using 
p-^-resummed instead of p-resummed flux. Lower panels: fractional difference between the numerical and analytical energy 
fluxes for the representative spin cases: q — —0.95, 0, 0.95. 



contrast, the agreement is substantially improved when 
we include the 5.5PN order nonspinning terms in the 
Taylor-expanded-truncated flux. The Taylor-expandcd- 
nontruncated flux agrees better with the numerical flux 
than the Taylor-expanded-truncated flux for retrograde 
orbits with q < 0, while its agreement is worse for pro- 
grade orbits with q > 0. For spin values q > 0.5, the 
agreement is especially bad, as the Taylor-nontruncatcd 
flux grows too fast when v > 0.4. We find that this differ- 
ence is mainly due to the large new spin term [l^ in the 
(3,3) mode, i.e. {-q"^ + 9tt q^ /2 + q89/5) v'^ in Z^s^^ 
(real part only). Without any resummation, the Taylor- 
expanded-truncated flux agrees well with the numerical 
flux for all spin values except for q = 0.95. The lower left 
panel shows that the fractional differences between the 
numerical and the Taylor-expanded-truncated fluxes are 
below 1% until v = 0.3, and are below 10% for q = 0.95 
until V = 0.55, and below 10% for all other spin values 
until the LSO. 

In the right panel of Fig. |9] we consider three analyt- 
ical flux models which use the fim, Pim (for q = 0.95 
only) and pj^, respectively. The fractional difference be- 
tween the numerical flux and /, p or p-'^-resummed fluxes 
is < 0.3% for all spin values when v < 0.3. Larger dif- 



ferences appear only when v > 0.3 for large and aligned 
spins, and the /-resummed flux performs especially bad 
when V > 0.4. In the case oi q = 0.95, we show the 
significant improvements achieved from the /-resummed 
to the p-resummed and eventually to the p-'^-rcsummcd 
fiux. The fractional difference with numerical fiux at the 
LSO is reduced from - 3.5 x 10^ to - 3 to 13%. The 
main reason for the bad performance of the /-resummed 
fiux is caused by the new spin term [l3| in the (3,3) 
mode, i.e. (— + 9 tt 17^/2 + g 89/5) v'' in ^33^^^ (real 
part only), as is in the case of the Taylor-expanded- 
nontruncated flux. As a matter of fact, we notice that if 
we did not include this new term computed in Rcf. [isj . 
and applied the /-resummation, or the p-rcsummation 
only to the nonspinning terms [l8l | j40l |. we would find a 
fiux not very different from the p-resummed fiux in the 
right panel of Fig. ^ In the p or p-^-resummation, this 
new term is suppressed by an order of magnitude, which 
leads to the improvements in their performance in mod- 
eling the numerical fiux. Specifically, this term becomes 
(g6187/330- g3) in /33, (g 5297/2970 + (7^3) in 

P33, and {-q 1073/1188 + q^ 2/3) in p^g. 

Finally, for large aligned spin q = 0.95 at the LSO, the 
p-^-resummed flux is closer to the numerical flux than 
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the Taylor-expandcd-truncatcd flux. Furthermore, we 
want to emphasize that the p-^-resummation improves the 
Taylor-expanded flux substantially over a large range of 
V and spin values. The differences between numerical 
and p-'^-resummed fluxes are smaller than those between 
the numerical and Taylor-expanded-truncated fluxes, by 
a factor of 3 - 5 at low velocities. Considering the large 
number of orbits an extreme mass-ratio binary spends in 
this range of velocities or frequencies, such an improve- 
ment is indeed si gnifi cant in correcting the orbital dy- 
namics (see Ref. [3l| for a quantitative analysis in the 
nonspinning case). 

IV. FACTORIZED MULTIPOLAR WAVEFORMS 
FOR GENERIC MASS-RATIO SPINNING, 
NON-PRECESSING BLACK HOLES 

In this section we extend the calculation of Sec. to 
generic mass-ratio spinning, non-precessing black-hole bi- 



naries. 

In Ref. [2l|, [l^, the non-spinning Taylor-expanded 
multipolar waveforms were computed through 3PN or- 
der. In Ref. (23j . spinning Taylor-expanded multipolar 
waveforms were computed through 1.5PN order. Using 
the definitions: 

M = nil + m2 (37a) 
Sm = mi — 7712 , (37b) 

XA^U^-^), (37d) 

2 \ 777^ 7772 / 

and restricting ourselves to circular, equatorial orbits, we 
obtain the following modes decomposed with respect to 
-2 spin- weighted spherical harmonics 
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The 1.5PN, 0.5PN and IPN order spin terms in the Ref. The 1.5PN-order (0.5PN-order) spin terms in 

modes /722, /121, /133, respectively, were obtained in the even (odd) parity modes are computed in Appendix 
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IfI The higher-order non-spinning PN terms can be found 
in Refs. [llllil. 

To compute the factorizcd multipolar waveforms for 
generic mass ratios we use Eq. ([T]) . For the source terms 

S^ff'' we employ the energy and angular momentum for 
circular, equatorial orbits computed from the effective- 
one-body Hamiltonian of Ref. [1^ (at the PN order 
at which we derive the factorized modes, the Taylor- 
expanded Hamiltonian of Ref. js^ coincides with the 
Hamiltonian of Ref. [s^ ) . More explicitly, when expand- 
ing the effective-one-body energy and angular momen- 
tum for circular, equatorial orbits through 1.5PN order, 
we find 
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Equations ([5^ . (^0]) are sufficient for computing the 
quantity fim in Eq. ([1]). In fact, similarly to the test- 
particle case analyzed in Sec. HIl the factor Tim in the 
generic mass-ratio case is not modified by spin effects. 
The factor Sim is not modified by spin effects either since 
there is no imaginary spin terms in Eqs. p8ap - (|38ip . 
The nonspinning 6em expressions for generic mass ratios 
are given in Eqs. (20)-(29) of Ref. p- Thus, inserting 
Eqs. ([39]) and gO]) in Eq. dTj), and using Eqs. (p8al) - ((38i| . 
we derive the even-parity f^m and p^m and odd-parity 
fern Pirn *° highest PN accuracy known to- 
day. We obtain 
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We may use E{r) instead of \L\ as the source term in 
the odd-parity modes. However, there is no difference 
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FIG. 10: Comparisons between the numerical, Taylor- 
expanded and p-resummed amplitudes of the dominant modes 
for an equal-mass equal-spin black hole binary as functions of 
the orbital velocity v. In the left panel, the component spins 
are Xi = X2 = 0.43655; in the right panel, the component 
spins are Xi = Xz = —0.43757. The numerical amplitudes 
were produced by the Caltech/Cornell/CITA collaboration. 



between ff"^ and /j^^, and correspondingly between pf^ 
and p/^, through PN orders where spin effects of binaries 
with generic mass ratio are known. 

In the nonspinning case, using IPN, 2PN and 3PN cor- 
rections, it was shown Q that the dependence of ptm on 
the mass-ratio v is mild. As a consequence, it was consid- 
ered meaningful to use test-particle results at PN orders 
where generic mass-ratio results are unknown. Since for 
each mode only the leading order generic mass-ratio spin 
terms arc known, it is not possible to carry out an ex- 
haustive study and understand how the spin terms in 
pim depend on v. As obtained in Appendix |F1 at lead- 
ing order, the 0.5PN spin terms in the odd-parity modes 
are proportional to v. Thus, they are zeros in the test- 
particle limit, but finite in the comparable mass case. 
Moreover, we find that the dependence on v of the 1.5PN 
spin terms in the even-parity modes is not that simple. 
Depending on the values of xs and XA, the relative differ- 
ence between h^i^^'^^^ {v = 0.25) and h^i^^'^^^ = 0) 
varies from zero to order of unity. Therefore, also the 
dependence of h^^^^^'^™ on v is not mild. 

Nevertheless, it is still reasonable to include the test- 
particle limit spin terms in fi„i and p^n such that at 
least part of the higher order spin effects are included, 
and to check the results against available numerical (ex- 
act) data. Specifically, we combine the test-particle and 
generic mass ratio results by replacing all the test-particle 
terms in /^m and p^m whose generic mass-ratio counter- 
parts are known with their generic expressions. 



Thus, in the generic mass-ratio, spinning case, we pro- 
pose to add to the fimS and pim^s derived in this section 
the test-particle limit terms derived in Sec. |lll In apply- 
ing this procedure we need to make a choice for the di- 
mensionless spin variable q appearing in the test-particle 
limit fim^s and pim^a. For a black- hole binary with com- 
ponent masses mi and m2 and spins xi and X2; we con- 
sider here two possibilities motivated by the choice of the 
deformed-Kerr-spin in the effective-one-body formalism. 
References [l^, [s^] used for the deformed-Kerr-spin 
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while Ref. [33[ used the following deformed-Kerr spin 
\S\ 1 



\/l-4i/XA + (l-2i^) xs- 



(44) 



Moreover, in the generic mass-ratio, spinning case, we 
also propose to use as effective sources in Eq. ([ij the 
Hamiltonian and angular momentum for quasi-circular 
orbits computed using the effective-one-body Hamiltoni- 
ans [Ulli. 

In Fig. 1101 we compare the amplitudes of the numeri- 
cal, the Taylor-expanded and the p-resummed modes for 
the five most dominant modes and for the two configu- 
rations of equal-mass, equal-spin black hole binaries of 
the Caltech-Cornell-CITA collaboration of Ref. 
We employ the effective sources built using the Hamilto- 
nian and angular momentum for quasi-circular orbits of 
Ref. .18, , j52| • The dimensionless spins in the two configu- 
rations are Xi = X2 = 0-43655 and Xi = X2 = —0.43757, 
respectively. The numerical amplitudes are derived from 
the numerical simulations published in Ref. [l^. Oscil- 
lations in the numerical amplitudes are due to numeri- 
cal artifacts in the simulations. For the (2,2) mode, the 
Taylor-expanded amplitudes agree quite well with the nu- 
merical amplitude, at least up to the frequency consid- 
ered. Thus, the improvement due to the p-resummation 
is marginal. For higher-order modes, there are large dif- 
ferences between numerical and Taylor-expanded ampli- 
tudes, and we find a substantial improvement when we 
adopt the p-resummation, except for the (3,3) mode in 
the spin aligned case (xi = X2 = 0.43655) whose numer- 
ical and Taylor-expanded amplitudes overlap, likely by 
coincidence. For the (2,2), (4,4) and (6,6) modes, the 
relative difference between numerical and p-resummed 
amplitudes is within 5% [4l|. For the (3,2) and (4,2) 
modes, the relative difference is between 10%-20%. We 
find that the results in Fig. [TO] depend weakly on the 
choice of q. In fact, using q = qo defined in Eq. (I43p and 
g = gs go/2 (when v = 0.25) defined in Eq. the 
relative amplitude difference is < 2% for the (2,2), (4,4), 
(4,2) and (6,6) modes, and ~ 5% for the (3,2) mode. 
Therefore, the uncertainty in the p-resummed amplitude 
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FIG. 11: Relative difference between p-resummed amplitudes 
of the dominant modes for an equal-mass, equal-spin black- 
hole binary when the test-particle spin is set to g = go and q = 
qs- The component spins of the binary are Xi = Xz = 0.95. 
The relative amplitude difference is plotted as a function of 
the orbital velocity v. 



due to the choice of q is less than half their systematic 
difference from the numerical results. 

Since we expect a stronger amplitude dependence on 
q in the case of larger spin magnitudes, we study the p- 
rcsummed amplitude dependence on the choice of q in 
Fig. [TTl where we show the difference between the am- 
plitudes \him{q = qo)\ and \herniq = qs)\ for an equal- 
mass, equal-spin black-hole binary with component spins 
Xi = X2 = 0.95. The relative amplitude differences are 

< 5% for the dominating (2,2) and (4,4) modes, and 

< 10% for the weaker (3,2) and (4,2) modes at t; < 0.45. 
For the (6,6) mode, since the test-particle spin terms in 
P66 a-re known only through 2PN order, i.e. only one more 
term is known beyond the generic mass ratio results, the 
amplitude dependence on q is entirely determined by this 
term and is somewhat stronger — reaching 30% at the 
LSO. 

Finally, we check the effect of the test-particle spin 
terms by comparing \himiq = qo)\ and \himiq = 0)| (i.e., 
removing the test-particle spin terms from the generic 
mass-ratio amplitudes in the latter) for this binary con- 
figuration. The difference, compared to Fig..[TTJ becomes 
larger by a factor of a few and reaches 10-25% for the 
(2,2), (4,4), (3,2) and (4,2) modes in the range of fre- 
quencies investigated in this paper. These terms may 
provide non-negligible corrections to the waveform and 
flux modeling. 



V. CONCLUSIONS 

In our study we employed the spin PN multipolar 
waveforms derived and decomposed with respect to the 
-2 spin- weighted spheroidal harmonics in Refs. [l3|, and 
transformed them in -2 spin-weighted spherical harmon- 
ics. We also took advantage of the new, recently com- 
puted (l3| . higher-order nonspinning and spin PN con- 
tributions in several subdominant modes. We also aug- 
mented our knowledge of the higher-order spin terms for 
generic mass-ratios, computing the generic expressions 
for the half and, one and half post-Newtonian contribu- 
tions to the odd-parity (current) and even-parity (odd) 
multipoles, respectively (see Appendix IF)) . 

Using the above results we extended the resummation 
method of factorized multipolar waveforms introduced in 
Ref. to spinning, non-precessing black-hole binaries. 
This factorized multipolar decomposition consists in a 
multiplicative decomposition of the him waveform into 
the product of several factors corresponding to various 
physical effects and the replacement of the factor ffm by 

its e-th root pim = {!lraYl^- 

In the case of a nonspinning test-particle orbiting a 
Kerr black hole in the equatorial plane, we found that 
the p-resummation is quite effective in reproducing the 
numerical multipolar amplitudes and energy flux up to 
q > 0.75 and v > 0.4. However, for larger values of g, we 
observed that the analytical pimivYs either have a slope 
larger than the numerical one, or they tend to grow as 
function of v instead of decreasing. This behavior can 
be cured by factoring out the lower-order PN terms in 
the pem, notably the 0.5PN, IPN and 1.5PN order terms. 
Being the lower-order PN terms negative (for q > 0), this 
procedure corresponds to factoring out the zeros of pim 
which turns out to capture the numerical (exact) zeros.. 

When applying the p-^-resummation, we found that 
the fractional amplitude difference between the numer- 
ical and analytical (2,2) mode at the LSO is 16% (33%), 
0.18% (0.32%) and 0.20% (0.85%) for q = 0.95,0,-0.95, 
respectively. Wc indicated in parenthesis the numbers 
when Taylor-expanded amplitudes arc employed. Thus, 
we found that for the (2,2) mode the improvement of the 
resummation is marginal. This might be due to the fact 
that the (2,2) mode is known at rather high PN order 
(5.5 PN for nonspinning terms and 4PN for spin terms). 
For the (2,1), (3,3) and (4,4) modes, for which less spin 
PN terms are known (see Table [T| the improvement due 
to the p-'^-resummation is even more striking. In fact, 
for those modes we obtained a fractional amplitude dif- 
ferences 2.4% (4.2), 0.2% (0.58%) and 0.0036% (0.15%), 
7.5% (2), 0.027% (0.55%) and 0.13% (0.2%), 16% (7.5), 
1.7% (28%) and 0.6% (5.8%), for q = 0.95,0,-0.95, re- 
spectively. For i > 5, the p-^-resummed amplitudes are 
certainly better than the Taylor-expanded amplitudes, 
but they differ from the numerical results quite substan- 
tially at high frequency. This is due to the fact that for 
those modes the spin effects arc known only up to 2.5PN 
order or lower. In summary, we found that the multipolar 
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amplitudes computed with the p-^-resummation are sys- 
tematicahy closer to the numerical (exact) results than 
Taylor-expanded ones over a large range of v and spin val- 
ues. The agreement can be further improved by including 
suitable adjustable parameters and calibrating them to 
the numerical results, as done in the non-spinning case 
in Ref. [U. 

Moreover, the numerical energy flux can also be suc- 
cessfully modeled by the p-^-resummation — for ex- 
ample we found that the fractional difference between 
the numerical and p-^-resummed flux is 13% (63%), 
0.70% (3.3%) and 0.48% (2.9%) for q = 0.95, 0, -0.95, re- 
spectively, where the numbers in parenthesis refer to the 
Taylor-expanded-truncated PN flux. For large aligned 
spins, the p-'^-resummed flux is much closer to the numer- 
ical flux at the LSO than the Taylor-expanded-truncated 
flux. Furthermore, we emphasize again that the - 
resummation improves the Taylor-expanded flux sub- 
stantially over a large range of v and spin values, and 
especially at low frequency where the majority of the 
signal-to-noise ratio of a binary accumulates. 

We have also extended the factorized resummation to 
generic mass-ratio, non-precessing, spinning black-hole 
binaries, and proposed, as in Ref. [6|, to augment the 
generic mass-ratio p^m with higher-order test-particle 
spin contributions. Unlike in the nonspinning case Q, 
in the spinning case only the leading-order generic mass- 
ratio spin terms are known. Using this limited informa- 
tion we found that the dependence on v of the spin terms 
is not necessarily mild. It depends on the mass ratio and 
the spin values. Nevertheless, we explored the possibil- 
ity of adding the spin contributions from the test-particle 
limit case to the generic mass-ratio amplitudes. 

When adding the test-particle limit contributions, we 
proposed to identify q with the Kcrr-deformed spin 
in the effective-one-body description. Using the two 
choices currently available in the literature, that is g = 
|So|/M2 [H il or g = jSj/M^ [H, ee found that the 
resummed amplitudes of the (2,2), (4,4), (4,2) and (6,6) 
modes agree with numerical simulation results |18| to 
within 2%, for equal-mass, equal-spin binaries with spins 
Ixil = 1x2 1 — 0.44. The (3,2) mode amplitude agrees 
with numerical results at 5% level. The relative differ- 
ence between the two choices of resummed amplitudes 
is less than half their difference from numerical results. 
When the spins are near extremal, e.g., Xi = X2 = 0.95, 
we found a mild, but non-negligible q dependence of the 
resummed amplitudes. Finally, when setting 9 = 0, that 
is removing the test-particle spin terms from the generic 
mass-ratio amplitudes, we obtain that the results vary by 
10-20% for the (2,2), (4,4), (3,2) and (4,2) modes in the 
range of frequencies investigated in this paper. 

The study carried out in this paper should be con- 
sidered as a first step in the modeling of extreme-mass- 
ratio inspirals and comparable mass black-hole binaries 
in presence of spins. We expect that in the extreme- 
mass-ratio inspiral case, the amplitude and flux agree- 
ment can be further improved by including in our p^^ a 



few adjustable parameters and calibrate them to the nu- 
merical data, as already done in Ref. j3l| for nonspinning 
extreme-mass-ratio inspirals. In the comparable-mass 
case, more detailed comparisons with accurate numerical- 
relativity simulations will allow us to nail down the choice 
of the spin parameter g, and allow us to carry out direct 
comparisons between the numerical and analytical p£m, 
thus helping in modeling the latter. 
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Appendix A: Taylor-expanded multipolar waveforms 

In order to compute the multipolar waveforms for a 
test particle around a Kerr black hole, we transform the 
Teukolsky equation into the frequency domain, and ex- 
pand it into the -2 spin-weighted spheroidal harmonics. 
The resulting equation is an ordinary differential equa- 
tion about the radial coordinate. This radial Teukolsky 
equation can be solved formally by using the Green func- 
tion. Since the Green function is represented by homo- 
geneous solutions of the radial Teukolsky equation, the 
central issue of this problem is to obtain the homogeneous 
solutions. There are two methods for obtaining them. 

In the first method, we transform the radial Teukol- 
sky equation into the Sasaki-Nakamura equation. In the 
Schwarzschild case, the homogeneous Sasaki-Nakamura 
equation becomes the homogeneous Regge- Wheeler equa- 
tion. We expand the homogeneous Sasaki-Nakamura or 
Regge- Wheeler equation in terms of e = GMu where lu 
is the angular frequency of the wave. In the case of cir- 
cular orbit, ui becomes ujq = mf2. (we revive the gravity 
constant G here). We look for the solution in power 
series in e. This is thus a kind of post-Minkowskian 
expansion. One difference between the ordinary post- 
Minkowskian approximation and this approximation is 
that we must impose correct boundary conditions at the 
horizon. Closed analytic representation of the solution at 
each order is necessary in order to obtain the asymptotic 
amplitudes which constitute the Green function. The 
lowest order solutions arc represented by spherical Bessel 
functions. The higher-order solutions can, in principle, 
be derived iteratively. However, it becomes more diffi- 
cult to perform this iteration and to derive the solution 
in closed analytic form at higher orders. The highest 
order computation so far was done in the Schwarzschild 
case by Tanaka et al. [ll| in which the closed analytic for- 
mulas for a homogeneous solution is obtained up to 0{e) 
for arbitrary £, and up to 0{e^) for £ = 2 and 3, and up 
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to ©(e^) for £ = 4. The formulas are explicitly given in a 
review paper j27| . Those computations are sufBcient for 
obtaining the energy flux through 5.5PN order. Since the 
formulas for the Zimujg 's are not given in the literature, 
we write them below. For each mode, we write the terms 
up to 0{v^^^^^^~'^'>) relative to the lowest-order term. 

Furthermore, in the Kerr case, so far the highest order 
computation was done by Tagoshi ct al. jl2| in which 
the closed analytic formulas for a homogeneous solution 
is obtained at 0(e) for arbitrary £ modes, and at O(e^) 
for £ = 2 and 3 modes. These computations are suflicient 
for obtaining the energy flux through 4PN order. 

Two of the authors have recently obtained the 0{e^) 
closed analytic formulas for i = 4 mode This order 
is necessary to derive the multipolar waveforms through 

3PN order beyond cl^^\ i.e. 4PN order beyond Cj^'^^ 
(see Table H]). More details of the computation and com- 
plete results are given elsewhere [isj . Here we show only 
the explicit formulas for Zf^wo defined in Eq. ^B)) . We 
write the spin-dependent 4PN-order Zimcjg in which each 
mode contains terms up to C'(t;^~^^~^-'~'^p) relative to the 
lowest-order term, (e^ is the parity of each mode). 
The second method to obtain the homogeneous Teukol- 



sky function is based on the Mano-Suzuki-Takasugi for- 
malism psj . In this formalism, the homogeneous solu- 
tions of the Teukolsky equation is represented with the 
series of hypergeometric functions and confluent hyper- 
geometric functions. The expansion coefficients of the 
two series solutions are the same, and they are closely 
related to the series expansion in power of e. Thus, if 
we compute this series up to higher order, we automat- 
ically obtain the higher order PN expansion formulas. 
Such computation was applied to the evaluation of the 
PN expansion of the black hole absorption effect in the 
Kerr case [IHl- This method was also applied to the en- 
ergy flux through 5.5PN order in the Schwarzschild case, 
confirming the results obtained with the above iteration 
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We apply this method to the Kerr case 
and obtain the 4PN-order multipolar waveforms which 
agree with the results obtained with the above iteration 
method. This method has also been recently applied 
to the computation of the 5.5PN-order multipolar wave- 
forms in the Schwarzschild case by Fujita and Iyer p^ . 
The non-spinning terms of expressions below agree with 
their results up to C(?;"-2(^-2)). 

We have 
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Appendix B: Expressions of the C^m's modes for 

4 < £ < 8 
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Appendix C: Expressions of the ftm's modes for £ > 4 
1. The odd-parity fi'^'s and even-parity /^m's 
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Appendix D: Expressions of the pcm's modes for £ > 4 
1. The odd-parity pf'^'s and even-parity ptm's 



31 



P55 = 1 - 



487 2 2(7 3 h 3353747\ , 241 g , . , , 

If -v^ + [- V* , (Dla) 

390 3 V 2 2129 400/ 195 ^ ' 

J 2908 2 2g . /2g2 16213384\ , ,^ , , 

- 1- 2275" -ll^' H^^I^^)^ ' 
2, _ £0833N 103, 

3 V 2 709 800/ 325 ' ^ ' 

2, /2^ _ 1187914 N 
15 V 5 15 526 875/ ' 



T 25 2 

pf, = l-g|.^-^.^+(^-V,^^).\ (Did) 



P51 



^ 319 2 2 (J 3 /q2 31877 \ 4 139 g 5 

1 if -V^ + [- H , (Die) 

390 3 V 2 304200/ 975 ' ^ ' 



53 2 2g 3 fq^ 1025435\ 4 l _i 185 2 2g 3 



42 3 V 2 659 736 ^ ' ''es 144 9 ' ^ ^ 

^ 43 2 2g 3 /q2 476887\ 4 . , 169 2 2g 3 . ^ , 



42 3 V 2 659 736/ 144 9 

37^2 2g 3 , fq^ 817991 ^4 ^ _^ 161 2 2g 3 



P62 = 1 2. + Pfii = 1 -v\ (Dlh) 

42 3 V 2 3 298 680 ^ ' ''ei 144 9 ' ^ ' 

1 151 2 2(7 3 , ^ 1072 2 , 
^ 127 2 2(7 3 , , 8878 2 

^-^i-m" -T" ' ^-=1-7497" ' ™ 

f^73 ^-^g 3 , P72 rj^gj , P71 119 3 ' ^ ^ 

. 1741 2 r . 3913 2 . 167 2 l . 725 2 .^..n 

P88 = 1 W , PsT = 1 V , P86 = 1 1' , Psr, = 1 , (Dll) 

1333 2 r , 3433 2 ^ 1 231 2 l . 3337 2 ^ 

P84 = 1 V , Ps-? = 1 , P82 = 1 V , Psi = 1 V . (Dim) 

F84 -^ggg , P83 3Q4Q ' P82 -^ggg , FBI 3 q^q V ^ 



2. The odd-parity p^^'s 
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Appendix E: Expressions of the (5fm's modes for 
4 < ^ < 7 
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Appendix F: Multipole moments for generic £ and m 

In Refs. [111 the authors have computed the even- 
and odd-parity IPN multipoles for generic £ and m. 
Those calculations were crucial in understanding the £- 
scaling of the fem^s, suggesting the introduction of the 
ogrnS functions. 



In this Appendix, we calculate the 0.5PN spin terms in 
the odd-parity multipoles h^^Ji ^'^'^ ^^'^ 1.5PN spin terms 
in the even-parity multipoles h^^Ji- J^^st for completeness 
we also reproduce the IPN nonspinning terms in the odd- 
parity multipoles h^i^, already computed in Ref. Q. 

Henceforth, we make use of the standard multi-index 
notation for tensors of arbitrary rank, which are dis- 
played as 



(Fl) 



where each index ii to Z£ runs from 1 to 3. We also 
employ the notation T<i> = STFL[Ti] to denote the 
symmetric trace-free projection over the indices ii to if. 
For example we have 



T. 
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1 



(F2) 



Repeated multi-indices imply summation over all corre- 
sponding indices, e.g. 



L*-* — ii i->...i,'-> 



(F3) 



Reference [2l| computed the expression of the full wave- 
form as an expansion in -2 spin-weighted spherical har- 
monics through the coefficients Ugm and Vim as follows 
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(F4) 
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where 



IGtt 
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Vir, 



327r£ 



(^ + 2) 



(2£ + i)!! V 2^{^ + l){^-l) 



(F5b) 

The radiative moments U l and Vl are the l^^ time deriva- 
tives of the multipole moments and respectively, 
as we neglect tail contributions for our purposes here. 
In terms of the vector f defined above Eqs. (jFlip . the 
quantity 3^£'" is defined as follows 



yrfL. 
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1. Odd-parity 0.5PN spin multipoles 



The odd-parity contributions to the waveforms are pro- 
vided by the expansion coefficients Vgm, which in turn 
are determined by the current multipole moments Jl ■ In 
the circular orbital case, the nonspinning IPN current- 
multipolc moment Jl is given by [6| 



tNS 



(F7) 



where fl is the orbital frequency, v ~ (A/fi)^/'^ and 
r 



-£n = -r-^^, A = iNXn, (F8) 
r X r 



and where 
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where C( coincides with Eq. and Xi^2 = 'mi,2lM. For 
circular orbits, we have 

n = (cos0orb,sin0orb,O) , (FlOa) 
A = (-sin(?!)orb,cos(?!)orb,0) , (FlOb) 
Ln = (0,0,1). (FlOc) 

In terms of the vector r — (sin6'cos </), sin6' sini^, cos 6'), 
the following expressions will prove very helpful below 
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The 0.5PN-order contribution to Jl that is linear in the 
spins is given by Ref. [sl] 

j|=(i±l)sTF,|^5-y^i|. (F12) 

To rewrite Eq. (|F12p in the ccnter-of-mass frame, we use 
Hi = X2 r and ?/2 = —Xi r, which leads to the following 
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(F13) 



where 



= XiXt'xi + i-Y-'X2Xt'x2, (F14) 

and we define xi = Si/rrii and X2 = S2/rn^. For non- 
precessing binaries, we have S(^) = ^n-, and hence 
we can write down the total IPN-order current multipole 
moment as 



W K2 n^"^ X'-^ 
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Next, in order to compute the radiative coefficient Vim, 



we first need 



dm 



It is therefore useful to 



rewrite all vectors appearing in in terms of r as follows 
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+V- 



L-l 



(F16) 



orb 



where the "orb" subscript is shorthand for evaluating the 
bracket at 6 — 7r/2, = t^orb- The purpose of this rewrit- 
ing is to allow us to eventually make use of Eq. (|F6p . 
together with the following identities 
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By substituting Eqs. (|F17|) into Eq. (|F16p . the current 
multipole moments become 



Jl = {vMr'^+^n)'&T¥L \ ^dgn 
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„ K2 or y(£+l)- r\ Contracting Eq. dFlSl with 3^^™ * then yields 
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(F18) 
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From the parity properties of associated Legendre poly- 
nomials, J^,„ is non-vanishing only if £ -|- m is odd. The 
next step consists of converting r^"*"^ into an expansion 
in V by means of Kepler's third law, 



and substituting it into Eq. (|F19p yields 
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Taking £ time derivatives and multiplying by the appropriate normalization factor finally gives 
327r£ r 
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The overall factor in front of the bracket in Eq. (|F23|) coincides with the Newtonian contribution as given by Eq. ([2]) , 
using Eqs. (jF4p and pb[) . Hence by definition [see Eq. ([T])], we find 
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Again, the IPN-order terms in the above equation were 
computed in Appendix A of Ref. @ . 

Quite interestingly, we find that in the nonspinning 
test-particle limit (m2 <C Tni,xi = = ami = 



q,X2 = 0), only the odd-parity mode £ = 2 contains 
the 0.5PN spin term, for all the other odd-parity modes 
the 0.5PN spin terms vanish. In fact, using Eqs. (jFQdp . 
(IF14P we find that if we set X2 = 0, the 0.5PN spin terms 
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reduces to 

f (1),0.5PN 



{£ + l)mlxi 



2[mi ni2 + (—1)^ Tn[ m2 



(F25) 



If ^ = 2, then h. 
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if ^ > 2, we have h 



(1),0.5PN _ 



(1).0.5PN 

Im 



= —3/2vq when i/ — ^ 0, while 
cc v qv and the latter goes 
to zero as ^ 0. The fact that the odd-parity modes 
with i > 2 vanish, is consistent with the -2 spin- weighted 
spherical C^m's computed in the main part of this paper. 
However, it is worth noticing that the odd-parity -2 spin- 
weighted spheroidal Z^^^s do contain 0.5PN spin terms. 

Moreover, for the case of finite symmetric mass-ratio 
I/, we find that the 0.5PN spin terms in Eq. (|F24I) co- 
incide with what was derived in PN theory [23- The 
^-dependence of the 0.5PN spin term in Eq. (|F24p varies 
depending on the binary mass ratio and the spin magni- 
tudes. For example we find that for maximally spinning 
and aligned black holes (xi = X2 = 1) if the masses are 
equal, the 0.5PN spin term in Eq. (jF24p scales as i, but 
if the masses are unequal, it generally doesn't scale as £. 
Finally, we derive the corresponding generic 0.5PN 



spin contributions to and p^^^- Since we know that 
there is no quadratic spin contribution at IPN order in 
/im' ^® need to introduce a IPN quadratic spin term in 
Pim- Thus, the spin portions read 



(1) 



The second contribution comes from the Newtonian mass 
multipole moments, in two different ways. First, since the 
coordinate transformation that takes us from a generic 
frame to the center-of-mass frame involves the spins at 
1.5PN order, the Newtonian mass multipole moments ac- 
quire a spin contribution when re-expressed in the center- 
of-mass frame. Second, when we use Kepler's law at 
1.5PN order to rewrite the orbital separation r as an ex- 
pansion in V = {Mfiy/^, spin terms are generated which 
contribute to the 1.5PN spinning waveform. In a general 
frame, the Newtonian mass multipole moments are given 
by 



= STFi 



m2 y2 



(F31) 



The coordinate transformation to the center-of-mass 
frame is given by (42j 
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em 
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2. Even-parity 1.5PN spin multipoles 



Vi = Xor ^ V X A . 



2/2 



AI 



■V X A . 



(F32a) 
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The 1.5PN spin contributions to the even-parity wave- 
form come from two distinct sources. The first is the 
1.5PN spin mass multipole moment /f , given by (in the 
center-of-mass frame, for non-prcccssing, circular orbits) 



where 
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Therefore in the center-of-mass frame, the Newtonian 
mass multipole moments read 



(F28) 



Making use of the following identity which is valid for 
circular orbits 



STF, 



(£-1) 



= STFl 



we can rewrite /f as follows 



(F29) 



If = Mvci,r^n<^> ^v"^ £ct-ir^-^n<^-^{v x A)*^> . 

(F33) 

For non-precessing, circular orbits, Eq. ()F33|) may be 
rewritten as 



op 

if = M'^v'^- — -/EmSTFi 



£+1 



(F30) 



if = Mvcir^n<^> 



£ci- 



Cl 



(F34) 
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Adding together both eontributions (|F30|) and (jF34p . 
contracting with and finaUy taking £ time derivatives 
as well as multiplying by the appropriate overall factor, 



we arrive at the following expression for the even-parity 
radiative moment 



16tt ii+l)ii + 2) 
(2£+l)!!V 2e{i-l) 



Miyci{-imnY/Y;„^{TT/2,M 



Ci 
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m 
1 



(F35) 



The overall factor in front of the bracket in Eq. (jF35[) 
coincides with the Newtonian contribution as given by 
Eq. ([2]), using Eqs. ((F4| and dSa]). Next we use Kepler's 
third law to replace the orbital separation r by the fol- 
lowing expansion in v. Aagain, we do not write the IPN 
order non-spinning contributions explicitly here to keep 
formulas short. 



r = Mt)"^ < 1+ 



';{xlxi+xlx2) + y{xi + X2) 



(F36) 

Substituting (|F36|) into (|F35p . we can finally isolate the 
1.5PN spin contribution to the even-parity waveform as 



- (0),1-5PN 



{Xl XI + Xl X2) +y{Xl+ X2) 



(^2X2-^1X1) 



2t 



(F37) 



In the non-spinning test particle limit, Eq. (|F37[) simply 
reduces to 



(0),1.5PN 



(0),1.5PN 
em 



(F39) 



h 



(0),1-5PN 



21 



(F38) 



thus, it scales as £. Finally, we derive the corresponding Therefore, the generic f^H and pfj^ expressions through 
generic 1.5PN spin contribution to /'"^ and pfj and they ^-^PN are given by the above equation combined with 
^gj^j thelPNnonspinningresuhgiveninEq. (A15) of Ref. @ 

(note that there is no 1.5PN nonspinning contribution to 

.(0),1.5PN _ r(0),1.5PN „(0) (0) 

Jem ~ "'iin ' Jem Pim' 
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